Coexistence of superconductivity and antiferromagnetism in self-doped bilayer t-t'-J 

model 



00 

o 
o 

(N 



I 

o 



> 

00 
00 

l> 

O 

oo 

o 



J. Y. Gan,i M. Mori,i T. K. Lee,^ and S. Maekawa^ 

'institute for Materials Research, Tohoku University, Sendai 980-8577 
^Institute of Physics, Academia Sinica, Nankang, Taipei, Taiwan 11529 
(Dated: July 17, 2008) 

A self-doped bilayer t-t'-J model of an electron- and a hole-doped planes is studied by the slave- 
boson mean-field theory. A hopping integral between the differently doped planes, which are gen- 
erated by a site potential, are renormalized by the electron-electron correlation. We find coexistent 
phases of antiferromagnetic (AFM) and superconducting orders, although the magnitudes of order 
parameters become more dissimilar in the bilayer away from half-filling. Fermi surfaces (FS's) with 
the AFM order show two pockets around the nodal and the anti-nodal regions. These results look 
like a composite of electron- and hole-doped FS's. In the nodal direction, the FS splitting is absent 
even in the bilayer system, since one band is flat due to the AFM order. 

PACS numbers: 74.72.Jt, 74.62.Dh, 74.20.-z, 79.60.-i 



I. INTRODUCTION 

High- Tc superconductors (HTSC) have one or more 
Cu02 planes in a conducting block, which is separated 
by charge-reservoir blocks. In HTSC with more than 
three Cu02 planes in a unit cell, there exist two inequiv- 
alent types of Cu02 planes; pyramidally-coordinated- 
outer planes (OP) and square-coordinated-inner planes 
(IP). The nuclear magnetic resonance (NMR) studies 
found that the hole density in OP is lager than that 
in IPi^i^iii^i^. An example of these is the five-layered 
HgBa2Ca4Cu50j,, in which the optimally doped OPs are 
superconducting (SC) with Tc = 108i^, while the three 
IPs have an antiferromagnetic (AFM) moment^iiii^. Al- 
though the SC planes are separated by the AFM ones, 
the Josephson coupling through the AFM planes stabi- 
lizes the superconductivity as a bulk^i. 

Another kind of multilayered HTSC is the four-layered 
Ba2Ca3Cu408(0:rFi-:,)2 (F0234)^. Especially for x = 
0, a nominal Cu valence is -1-2 on the canonical chemical 
formula. Thus, this material is expected to be a Mott 
insulator, although the superconductivity with Tc=60K 
takes place^ii^. This compound, F0234, has four Cu02 
planes, among which two OP's have apical F atoms, 
while the two IP's do not. Angle-resolved-photoemission- 
spectroscopy (ARPES) experiments observed two Fermi 
surfaces (FS's), whose volumes in the first Brillouin zone 
correspond to electron- and hole-doped FS'si2,iii. This 
would be the first self-doped high Tc superconductor with 
an electron- and a hole-doped Cu02 planes in the same 
crystal. It is also found that the superconducting (SC) 
gap on the electron-doped FS is twice as large as that on 
the hole-doped oneiS. 

On the other hand, it is known that doped holes make 
a FS around the nodal regioni^ii^iiilii^, while doped elec- 
trons create pockets around the anti-nodal regionsi^. 
Theoretical studies by the variational Monte Carlo 
mcthodiiii^and the exact diagonalization metho d^^i^° 
elucidate that the asymmetry between hole- and electron- 



doped cuprates results from second neighbor hopping (t') 
and third neighbor one {t") in the Cu02 plane. Here, the 
question arises; what is the ground state of the self-doped 
bilayer cuprates, where one plane is electron-doped and 
the other is hole-doped, and how are the FS's and their 
asymmetryl 

To answer these questions, the self-doped t-t'-J model 
is examined by the slave-boson mean-field theory. The 
two different types of planes are connected by an inter- 
layer hopping renormalized by electron-electron correla- 
tion. A site potential making the charge imbalance be- 
tween two planes is included. Note that the hopping of a 
single spin between a holon- and a doublon-sites picks up 
extra minus sign as compared to that between a holon- 
and a single-occupied sitesi^. This doublon effect leads 
to a spin singlet states between the two planes. 

In an undoped bilayer system, both planes have same 
amount of carriers due to the self-doping, although one 
type of carrier is hole and the other is electron. Our re- 
sults show that in both electron- and hole-doped planes, 
AFM and SC coexist. Other authors studied the un- 
doped case with no interlayer hopping^ii^. We examine 
the doped case with holes as well. In the doped case, 
the numbers of carriers in each plane becomes imbal- 
anced, i.e. doublon density decreases and holon density 
increases with hole doping. As a result, the magnitudes 
of order parameters become more dissimilar compared to 
the undoped case. Two FS's in the self-doped bilayer look 
like a composite of hole- and electron-doped cuprates. 
However, we cannot find the FS splitting in the nodal 
direction, since one band becomes very flat due to the 
AFM orders. 

This paper is organized as follows. In Sec. |TT1 we in- 
troduce the bilayer t-t'-J model with an interlayer hop- 
ping and a site potential, and present the slave-boson 
mean-field scheme. In Sec. IIIIl we discuss self-consistent 
mean-field solutions for both undoped and doped cases 
of self-doped bilayer system. Coexistent phase of SC and 
AFM orders are discussed from the viewpoint of dop- 
ing and charge imbalance. FS and dispersion relation of 
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spinons are shown in the AFM ordered phase. In Sec.lIVI 
we will give summary and discussion. 



II. MODEL AND METHOD 

The model we apply to study the self-doped bilayer 
system is written as: 



H = H\\+ Hw + 
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nf^),, where r] = x,y indicates 



the nearest-neighbor sites. Although the magnetic or- 
der in real materials may be quite complicated, we only 
consider the commensurate antiferromagnetic orders for 
simplicity. As the interlayer hopping may induce a weak 
AFM correlation between the two planes, the staggered 
AFM order has a sign difference between the two planes 
in our definition. All parameters are assumed to be real 
and the SC pairing symmetry is d-wave. 

The Hamiltonian (1) based on the above treatment is 
decoupled in the momentum space as follows: 
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where cf^ {cta^) is the electron annihilation (creation) 
operator with spin a at site i in the /-th plane. The 
electron number in each plane is denoted by n^P = 
X^o- CiCT^c[^\ and the averaged electron density is defined 
as, n = (n'-'^' + 'nP'^)/2. The signs, {ij) and {ij), run 
over nearest- and next-nearest neighbor sites, respec- 
tively. The chemical potential ^ and the site potential W 
control the charge imbalance. Below, we take J/t = f/3 
and t'/t = -0.4. 

The interlayer hopping in Eq. ^ has the dis- 
persion relation in the momentum space, £j_^fc = 

(i_L/4) (cosfc X COS ky^ , where tj_ is the amplitude with- 
out renormalization22i2ii^i2^i^. 

We treat Hamiltonian ^ in the slave-boson mean- 
field theory. The electron operator is represented as, 
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with /if ■* and being the 



bosonic holon and doublon operators, respectively^^iS^. 
The fermionic spinon operator is denoted by ff^ . In 
the self-doped case, we assume that one plane is hole- 
doped and the other is electron-doped. For the hole- 
doped plane, as there is no doublon, the electron op- 



erator can be expressed as, c. 



(0 



ffJhP\ with the 



Ea£^^fla = 1, while for the 



constraint, h^'^ 
electron-doped plane, as there is no holon, the electron 
operator can be expressed as, = o'f^l^^df\ with the 

constraint, df^^df^ + J2a fia'' fia — ^- Since we are 
interested in the electronic states at low temperatures, 
the boson condensation is assumed in each plane, i.e.. 



(Ot\ 



and (df^> = (df^^) 



where (5^'^ and S'"J^' are the holon and doublon densities. 

To decouple the Hamiltonian, we introduce the order 
parameters in the electron- and the hole-doped planes 



as, A^'' 
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where 7^ = 2(cosfca; + cosfcj,), rjk ~ 2(cosA:j; — cos fey), 
~ 4 cos kx cos ky , and k runs over the magnetic Bril- 
louin zone with \kx i ky\ < tt. Q = (7r,7r) is the mag- 
ne tic vect or and A*" is the total number of lattice sites. 

\fs^p's^ is the rcnormalization factor of t±—. Here, we 
assumed the 1=1 (2) is the hole (electron) doped plane. 

We note that as seen in the fourth term in Eq. ([5]), 
the interlayer hopping in self-doped bilayer system may 
induce an interlayer singlet-paring, which can be defined 
as,A,.(/W/(^)-/W/if). 

The momentum dependence of dispersions is given by 
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where 
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'EafaU and 
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J2a are spinon 



densities in plane 1 and 2, respectively. For the given 
total electron number n and site potential W, the mean- 
field parameters A^'^ x^^^ and m^^\ the charge density 
in each plane Sj^'' i^li ^) and the chemical potential fi are 
self-consistcntly determined in numerical calculations. 



III. RESULTS AND DISCUSSION 

A. Self-doped bilayer at half-filling (71=1) 

First, we focus our study on the undoped case, i.e., 
n = 1. In this case, the holon density in hole-doped 
plane is equal to the doublon density in electron-doped 



plane, i.e., S 
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Figure [T] shows the results of the 
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d-wave pairing amplitude (A^'^), the uniform bond order 
parameter (x'''), the AFM order parameter (m'^'^), the 
site potential (W) and the interlayer singlet pairing am- 
plitude (Ap) as functions of the holon (doublon) density 

S^^^ (^d^' ) fo'^ various values of the interlayer hopping pa- 
rameter {t±). A('\ x^'), and m*^'-* depend very weakly on 



t 



, particularly for small W (small S/^ 



(1) 



and^i')) 



Ap in- 



creases with t±. When t± = 0, Ap = 0. Fig.[T](a) shows 

the relation between W and Sj^^ (^i^')- When W = 0, 
there is no charge imbalance between the two planes, i.e., 
0. In this case, both planes are at half- 



^(1) = ^(2) 



filling, with m*-'^ 



1 and A(') = x^'^ = 0. The ground 



state is an AFM insulator. When (5^^^ — 



0, Ap be- 



comes zero, that is, the planes are decoupled regardless of 
the interlayer hopping and the planes are coupled only for 
finite J^^"* When W increases, the charge densi- 

ties 5^f^^ and increase from zero. The staggered AFM 
magnetization decreases with 5^^^^ {^^d^)i while the d-wave 
paring amplitude (A*^')) and the uniform bond order pa- 
rameter (x*'''') both increase. rrS^^ and to'^' are almost 



the same in the region < (5 



(1) r(2) 



< 0.15. and then 



77i'^^^ decreases faster than ra'^'^^ and vanish at around 
5\l^ = ~ 0.2. It is seen that for < ^ 0.2, 

both electron- and hole-doped planes are the coexistent 
state of AFM and SC. When AFM order vanishes, both 
planes are superconducting. 



(1) 



Here, we mention yet another solution, where 
S^^'> ^ even for W = 0, and m'^) > m^i) ~ 0. This 
may be a possible phase separation, where electrons are 
spontaneously transferred from the hole-doped plane to 
the electron-doped one to gain an energy of magnetic 
exchange interaction without the site potential. Since 
this solution is found in a limited (unphysical) parameter 
region, we do not discuss below. 
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FIG. 1: The order parameters in the undoped case with t' /t — 
—0.4. The antiferromagnetic (AFM) and the superconducting 
(SC) orders coexist in both planes. The following properties 
in electron- and hole-doped planes are plotted as functions 
of holon (doublon) density; (a) site potential, W, (b) d-wave 

(d) 



pairing amplitude, A''\ (c) the uniform bond order, 



AFM order parameter, m^'\ and (e) interlayer singlet-paring, 
Ap. Three lines in each figure correspond to different value 
of interlayer hopping given by, t±/t=0.0 (square), t±/t=0.5 
(circle), t±/t=0.0 (triangle). 



B. Hole-doped case {n <1) 

Next we investigate the doped case in the self-doped 
bilayer system. In the doped case, i.e., n 1, the holon 

density (Sj^^' is not necessarily equal to doublon density 
S'-^K Figure m shows the results of S^^ ,5f\ A^), x^'\ 
and TO*-'^ as functions of the total electron density n for 
a given W/t = 0.05. From Fig. [2] (a), we see that for 

small doping of holes, 0.97 < n < \, both (5^^^' and (5^^' 
increase. This means that holes first go into hole-doped 
plane, while some electrons arc transferred from hole- 
doped plane to electron-doped plane. The kinetic energy 
gains in this case. Upon further doping of holes into the 
system, holes go into both planes, and 5^^^ increases while 
6\ decreases. Due to the change of charge density, the 
staggered AFM magnetization in hole-doped plane (m'^' ) 
decreases while m'-^-' increases. For 0.86 < n < 1, both 
electron- and hole-doped planes are the coexistent state 



of AFM and SC. For 0.86 < 7i < 0.97, the AFM order in 
hole-doped plane is small, and m^^^ decreases fast with 
When = 0, nS^^ becomes zero in this region. 
At a critical point n ^ 0.86, doublon vanishes [S'^p = 0), 
and electron-doped plane goes into AFM insulator phase; 
simultaneously rrS^'^ becomes zero, and hole-doped plane 
goes into the superconducting phase. Above the critical 
point [n < 0.86), both planes are hole-doped. 

So far we have presented results for both undoped and 
doped cases. Now we discuss the phase diagram. Fig. [3] 
shows the phase diagram in the W-n plane for t±/t = 0.5. 
The phase diagram is divided into three parts. For small 
doping {n close to 1) and large site potential W, the 
charge imbalance is large and both planes are SC; for 
small W and small n, both planes are hole-doped; for 
intermediate W and n, both electron- and hole-doped 
planes are the coexistent state of AFM and SC. In the 
undoped case {n — 1), with increasing W, S'j^ = S'^p 
increases and there is a transition from coexistent state 
of AFM and SC to SC in both planes. In the doped case 
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FIG. 2: The order parameters in the doped case with 
t' /t = —0.4 and W/t = 0.05. The foUowing properties in 
electron- and hole-doped planes are plotted as functions of 
total electron number n; (a) holon (doublon) density, (b) d- 
wave pairing amplitude, A*-'', (c) the uniform bond order, 
and (d) AFM order parameter, m''\ Three lines in each 
figure correspond to different value of interlayer hopping given 
by, t±/t=0.0 (square), t±/t=0.5 (circle), t±/t=0.0 (triangle). 



with n ~ 0.85, both planes arc holc-dopcd for small W. 
When W increases above a critical value of Wc, electrons 
move from hole-doped plane to electron-doped plane and 
both planes are the AFM and SC coexistent state due to 
self-doping. 



C. Spinon Fermi surfaces and dispersion relation 

As for the asymmetry between the hole- and the 
electron-doped cuprates, one of distinguished observa- 
tions is the FS pocket, which is located around the 
nodal region in hole-doped cuprate a^^'^'^i^^'^^ and the 
anti-nodal region in electron-doped oncsiS. It is found 
that this asymmetry originates from the different signs 
of t' and t"2U^^l^. On the other hand, the multilay- 
ered cuprates doped with holes show the interlayer split- 
tings of FS^i^i25^i27^^^_ ^Yic nodal direction, the 
splitting means the charge imbalance between IP and OP, 
while those around the anti-nodal regions are ascribed to 
a magnitude of interlayer hopping renormalized by the 
charge imbalance^. Interesting is that the two asym- 
metric planes are combined by the interlayer hopping in 
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FIG. 3: The phase diagram in the W-n plane for t±/t = 
0.5. Each abbreviation means the foUowings: AFM, the an- 
tiferromagnetic phase: SC, d-wave superconducting phase: 
AFM-I-SC, the coexistent phase of AFM and SC. For large 
n, both planes are hole doped. Different lines correspond to 
different values of t' /t as, -0.2, -0,3 and -0.4 from bottom to 
top. 



the self-doped bilayer system. 

In Fig. m the spinon FS's and dispersion relations of 
self-doped bilayer systems arc plotted for some doping 
rates in the AFM phase. Details of parameters are in- 
cluded in the caption of Fig. U) Near the half-filling 

given by, n=0.98, =0.130, and =0.090, two FS 
pockets appear in the nodal and the antinodal regions 
shown in Fig. [4] (a). It looks like a composite of hole- and 
electron-doped cuprates. As shown in Fig.[3](e), four dis- 
persion relations are separated each another due to the 
large AFM moments, m(i)=0.47 and to(2'=0.68. With 
increasing hole density in the bilayer as n=0.95 and 0.9, 
the hole-doped like FS becomes larger as shown in Fig. 3] 
(b) and (c), and the AFM moment in the hole-doped 
plane markedly becomes small, to^^-'=0.05 and 0.01. As 
a result, the separation of four bands close on each an- 
other as shown in Figs. 2] (f) and (g). On the other hand, 
since the AFM moment in the electron-doped plane still 
large, to'^'=0.73 and 0.95, two among four bands become 
quite flat. Finally, for 71=0.85, both plane becomes hole- 
doped. Interesting is that we cannot find the interlayer 
sphtting in the nodal direction as shown in Fig. [5] (d), 
although it is found in the normal metallic phases^^i^. 
The missing of FS splitting is caused by the AFM mo- 
ment in the 2nd plane, which make a band flat as shown 
inFig.SKh). 
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scribed by the t-t'-J model, and the interlayer hopping 
and a site potential are included. In an undoped bilayer 
system, both planes have same amount of carriers due 
to the self-doping, although one type of carrier is hole 
and the other is electron. Our results show that in both 
electron- and hole-doped planes, AFM and SC coexist. 
In the doped cases with holes, the numbers of carriers 
in each plane becomes imbalanced, i.e. doublon density 
decreases and holon density increases with hole doping. 
The magnitudes of order parameters become more dis- 
similar compared to the undoped case. At some critical 
doping of holes, the doublon disappears and the electron 
layer becomes an insulator. This effect might be useful 
for a p-n junction^^ made of the electron- and hole- doped 
layers if we could control the doping near this value. Two 
FS's in the self-doped bilayer look like a composite of 
hole- and electron-doped cupratcs. However, we cannot 
find the FS splitting in the nodal direction even in the 
bilayer system, since one band becomes very flat due to 
the AFM orders. 

In the ARPES experiment on F0234, two FS's sur- 
rounding the (7r,7r) point and the FS's splitting in the 
nodal direction are observed. The doping rate in each 
plane may be optimal or overdoped. In such a case, the 
FS should enclose the (7r,7r) point like as Bi-compounds, 
not like as a pocket. On the other hand, the NMR study 
observed the magnetic moment, which could not exist 
around the optimum doping region. Although our model 
is a bilayered system, it involves essential points of mul- 
tilayered cuprates. In addition to the bilayered system, 
the four-layered t-t'-J model was examined to find the 
FS splitting in the nodal region. However, we could not 
find it in the self-doped four-layered system, while it was 
found in the four hole-doped layers. The contradiction 
between experiment and theory remains to be resolved 
in the future. 



FIG. 4: FS's and dispersion relations of self-doped bi- 
layer system; (a) and (e) n=0.98, (5{^^'=0.130, (jf '=0.090; (b) 
and (f) n=0.95, (5^^^=0.180, (jf^ =0.080; (c) and (g) n=0.90, 
5^^'=0.225, (5*,^'=0.025; (d) and (h) n=0.85, 5j^^'=0.25, 
5^^^=0.05. In these figures, W/f=0.05 and fx/t=0.5 are fixed, 
and SC order parameters are imposed to be zero. In (d) and 
(h), both two planes are in hole-doped regions. 

IV. SUMMARY AND DISCUSSION 

We have studied the bilayer self-doped cuprates by us- 
ing the slave-boson mean-field theory. Each plane is de- 
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